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The B meson decay modes B — > Dpp and B — > D*pp have been studied using 29.4 fb _1 of data 
collected with the Belle detector at KEKB. The B° — > D°pp and B° — > D*°pp decays have been 
observed for the first time with branching fractions B(B° — > D°pp) — (1.18 ± 0.15 ± 0.16) x 10~ 4 
and B(B° -> D*°pp) = (1.20±gJ| ± 0.21) x 10~ 4 . No signal has been found for the B+ -> D+pp 
and B + — » D* + pp decay modes, and the corresponding upper limits at 90% C.L. are presented. 
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To date, information on B decays with baryons in the 
final states is rather scarce. A recent search for the two- 
body baryonic B decays by Belle showed that their rel- 
ative probabilities are rather small: for the decay modes 
B° — » pp, AA, and B + — > pA upper limits of (1-2) x 1CT 6 
were obtained (!]]. At the same time, measurements of 
the B° — » D*~ppir + and B° — > D*~pn decay branch- 
ing fractions by CLEO (2| and the observation of the 
_B + — > ppK + decay by Belle S indicate the dominance 
of multibody final states in decays of B mesons into 
baryons. Moreover, the observation of the B° — > Z?°7r°, 
5° -> L>°r/ and B° -> L>°w decays by Belle (§ and 
CLEO || with branching fractions considerably higher 
than expected indicates that color-suppressed B decays 
with baryons in the final state are likely to be sizeable. 
This motivated a search for the color-suppressed decays 
B — > D^pp. The inclusion of charge conjugate modes is 
implicit throughout this report. 

We use a data sample collected with the Belle detector 
at the KEKB asymmetric energy e + e~~ collider f|. It 
consists of 29.4 fb _1 taken at the T(4S I ) resonance corre- 
sponding to N B g = 31.9 x 10 6 produced BB pairs, and 
2.3 fb _1 taken 60 MeV below the T(45) resonance to per- 
form systematic studies of the e + e~ — > qq background. 

The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a three-layer silicon ver- 
tex detector (SVD), a 50-layer central drift chamber 
(CDC) for charged particle tracking and specific ioniza- 
tion measurement {dE/dx), an array of aerogel thresh- 
old Cerenkov counters (ACC), time-of- flight scintillation 
counters (TOF), and a CsI(Tl) electromagnetic calorime- 
ter (ECL) located in the magnetic volume. The magnetic 
field is returned via an iron yoke that is instrumented to 
detect muons and Kl mesons (KLM). 

Charged tracks are selected with requirements based 
on the average hit residual and impact parameter rela- 
tive to the interaction point. We also require that the 
transverse momentum of the tracks be greater than 0.1 
GeV/c to reduce the low momentum combinatorial back- 
ground. 

For particle identification (PID), the combined in- 
formation from CDC, TOF and ACC subsystems is 
used. Protons and antiprotons are selected with a set 
of PID criteria that has an efficiency of 98% and a kaon 
misidentification probability of 15%. Selection criteria 
for charged kaons provide an efficiency of 88%, a pion 
misidentification probability of 8%, and negligible con- 
tamination from protons. All tracks positively identified 
as electrons are rejected. 

A pair of calorimeter showers with an invariant mass 
within 15 MeV/c 2 of the nominal ir° mass is considered 
as a 7T° candidate. An energy of at least 50 MeV and a 
photon-like shape are required for each shower. 

We reconstruct D mesons in the following decay chan- 
nels: D° -> K-TT+, D° -> K-ii+iT+ir-, D" -> K-n+n" 
and D + — » K~it + tt + . We select D candidates using 



cuts around the central values of the M(D) distribu- 
tions that correspond to 95% efficiency. For the ir° from 
the D° — > K~ir + Tr° decay, we require that the 7r° mo- 
mentum in the T(45) center-of-mass (CM) frame be 
greater than 0.2 GeV/c in order to reduce combinato- 
rial background. D* mesons are reconstructed in the 
D*° -> D°ir° and D*+ -> D°ir+ decay modes. Since the 
pions from D* — > Dir decays are slow, we relax these cuts 
and impose an energy threshold for tt° photons of 30 MeV 
and a transverse momentum threshold of 50 MeV/c. 
The mass difference between D* and D candidates is re- 
quired to be within 4 MeV from the expected value for 
D*° and 2.5 MeV for D*+ (~ 3er in both cases). 

We combine I)W candidates with pp pairs to form B 
mesons. Candidate events are identified by their CM 
energy difference, AE = (£\ Ei) — Eb, and the beam 
constrained mass, M bc = y/E'fi — (X^Pi) 2 ) where Eb is 
the beam energy and pi and Ei are the momenta and 
energies of the decay products of the B meson in the 
CM frame. We select events with M bc > 5.20 GeV/c 2 
and |A.E| < 0.2 GeV, and define a B signal region of 
5.272 GeV/c 2 < M bc < 5.288 GeV/c 2 and \AE\ < 
0.020 GeV. In the cases when there is more than one 
candidate in an event, the B — > Dpp (or B —> D*pp) 
candidate with the D mass (or D* — D mass difference) 
closest to the world average is chosen. We use Monte 
Carlo (MC) simulation with a three-body phase space 
distribution for the B — > D^pp decays to model the 
response of the detector and determine the efficiency [§J . 

To suppress the large combinatorial background that 
is dominated by the two-jet-like e + e~ — > qq continuum 
process, variables that characterize the event topology 
are used. We require |cos6* t hr| < 0.80, where # t hr is the 
angle between the thrust axis of the B candidate and 
that of the rest of the event. This cut eliminates 77% of 
the continuum background and retains 78% of the signal 
events. We also define a Fisher discriminant, J-, which 
includes: the production angle of the B candidate, the 
angle of the B candidate thrust axis with respect to the 
beam axis, and nine parameters that characterize the mo- 
mentum flow in the event relative to the B candidate 
thrust axis in the CM frame @. We impose a require- 
ment on T that rejects 28% of the remaining continuum 
background and retains 95% of the signal. 

The AE distributions were fitted with a Gaussian for 
the signal and a linear function for the background. The 
Gaussian mean value and width were fixed from the MC 
simulation of the signal events. In the fit to the AE 
distribution, the region AE < —0.1 GeV is excluded to 
avoid contributions from other B decays, such as B — > 
D^npp. For the calculation of branching fractions, we 
use the signal yields determined from the fit to the AE 
distribution. This minimizes a possible bias from other 
B meson decays, which tend to peak in M bc but not in 
AE. 

The AE distributions for the B° -> D°pp and B° -> 
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FIG. 1: The AE distributions for the B° -> D°pp candidates: 
(a) D° -> K~ir + , (b) D° -» K~tt + -k + tt- and (c) D° -> 
K~iy + ty°. The points with errors are experimental data, the 
hatched histograms are D° mass sidebands and the curves are 
fit results. 



D*°pp decays are shown in Figs. |l| and 0. The fit re- 
sults are presented in Table B Statistically significant 
signals are observed for the B° — > D°pp decay mode 
in all three decay channels of the D° meson. The cor- 
responding branching fractions are in good agreement 
with each other. For the final result we use a simul- 
taneous fit to the three D° decay channels. The AE 
distributions for each mode were fitted together by a 
sum of a signal Gaussian and linear background func- 
tion taking into account the corresponding detection ef- 
ficiencies and D° meson branching fractions. The nor- 
malization of the background was allowed to float while 
the signal yields were required to satisfy the constraint: 
Ni = N BS ■ B{B° -> D°pp) ■ B{D° -► X,) ■ e h where 
the branching fraction B(B° — ► D°pp) is a fit parameter; 
N BB is the number of BB pairs B(D° -> X t ) are 
the D° meson branching fractions to the final states Xi 
and Ci are the corresponding efficiencies. 

The signals in the B° —> D*°pp decay mode are less 
prominent but, when combined, have a 5.6a statistical 
significance. As a cross-check, we confirm that the dis- 
tribution in the -0.2 GeV< AE < -0.15 GeV region of 
the B° — > D°pp mode is consistent with background from 
B° — > D*°pp with this measured branching fraction. 

The B + — > D + pp and B + — > D* + pp decays are doubly 
CKM suppressed and, thus, are expected to have much 
smaller branching fractions. This is confirmed by the 
analysis of the corresponding distributions: we do not ob- 



FIG. 2: The AE distributions for the B° -> D*°pp candi- 
dates: (a) D° -> R--K+, (b) D° -» K--K+TT+-K- and (c) 
D° -» if" 7r + 7r . The points with errors are experimental 
data and the curves are fit results. 



serve any signal for the B + — > D + pp and B + — > D* + pp 
decays and present for them upper limits at 90% confi- 
dence level (C.L.) The Feldman- Cousins procedure |ll[] 
was used to calculate upper limits except for the simul- 
taneous fit, where the maximum likelihood method was 
applied. In this case the upper limit N was calculated 
from the relation J Q C{n)dn = 0.9 J °° jC(n)dn, where 
C(n) is the maximum likelihood with the signal yield at 
n. The systematic uncertainties were taken into account 
in these calculations. 

Figure |^ (a) shows the Dalitz plot for the B° — > D°pp 
candidates from the B signal region. For comparison, 
also shown in Fig. || (a) is the same distribution for MC 
B° — > D°pp signal events generated according to phase 
space. It is worth noting that apart from a threshold en- 
hancement in the invariant mass of pp (and possibly also 
D°p), the main part of the signal is distributed accord- 

* v pp 



ing to phase space. The Dalitz plot for the B° — > D 
channel (not shown) with a smaller statistics also reveals 
a similar tendency. 

Since the pp invariant mass distribution of the observed 
signal is not completely described by the phase space dis- 
tribution and the detection efficiency can be non-uniform 
over the Dalitz plot, some systematic uncertainty in the 
efficiency calculations may occur. To study the model 
dependence of the branching fractions, we fit the AE 
distribution for B° — > D°pp candidates in six bins of pp 
invariant mass and calculate the partial branching frac- 
tion separately for each bin. The results are presented 
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TABLE I: Branching fractions and 90% C.L. upper limits for B — * D^*'pp decays. 



Mode 








AA yield 


Af 6c yield 


Efficiency, % 


B, 10" 4 


Significance 


B° 


_> 








33.6±^ 


34.5l 6 6 | 


23.56 ±0.49 


1.17l» J 3 ± 0-14 


8.9(7 


B° 


-» 


D°pp, £>° -» A'-TT+Tr+Tr" 






24.2±|;? 


1<L7±\\ 


7.11 ±0.21 


1.42±S;|I ± 0.22 


5.6<r 


B° 


-> 


D°pp, £>° -» ^-tt+tt 






34.2±*j 


36.5t?;4 


7.28 ±0.28 


1.06±^12 ± 0.15 


5.1(7 


B° 


-> 


D°pp, simultaneous fit 






— 


— 


— 


1.18 ±0.15 ±0.16 


12(7 


B° 




£>*°pp, £>*° -» D°->K-tt 


+ 


5.0t 2 2 . 8 2 


e.i±i:i 


8.11 ±0.30 


0.8llg; 4 6 ±0.13 


2.9cr 


B° 


-» 


D*°pp, D*° D° -+ if-7r 


+ 7T+7T- 


3.5± I; 4 




1.96 ±0.15 


1.2llg;^ ±0.23 


2.6cr 


B° 




£>*°pp, D*° - DV\ D° - 




+ 7T 


I0.8l|;2 


13.6±|;| 


2.38 ±0.16 


1.65t°;g±0.30 


4.2cr 


B° 




D*°pp, simultaneous fit 












1.20±^;|| ±0.21 


5.6cr 


B+ 


> 


D + pp, D + -> K-TT+TY+ 






< 5.2 


< 5.1 


14.28 ±0.38 


< 0.15 90% C.L. 




B+ 


> 


D*+pp, D*+ -> D°tv+, D° 


-f a~ 


7T+ 


< 2.2 


< 2.3 


9.55 ±0.31 


< 0.34 90% C.L. 




B+ 




D*+pp, D*+ -> D°tv+, D° 


-> A'" 


7r + 7T + 7T _ 


< 1.8 


< 2.4 


2.54 ±0.16 


< 0.53 90% C.L. 




B+ 




D* + pp, D*+ -> D°tv+, D° 


-> K~ 


7T+7r° 


< 4.8 


< 6.2 


3.18 ±0.18 


< 0.61 90% C.L. 




B+ 




D* + pp, simultaneous fit 












< 0.15 90% C.L. 





TABLE II: Branching fraction for the B° — > D°pp in bins of 
the pp invariant mass. 



M(pp), GeV 


AA yield 


Efficiency, % 


B, 10" s 


< 2.13 


29.4+_^ 6 9 


10.01 ±0.44 


3.65±^3 ±0.51 


2.13 - 2.38 




9.01 ±0.32 


1.60lg;^ ± 0.22 


2.38 - 2.63 


13.4± 4 4 * 


8.15 ±0.28 


2.04±°:6i ± 0-29 


2.63 - 2.88 


14.51" 


9.57 ±0.30 


l.SSlo.S ±0.26 


2.88-3.13 


15.4±|;i 


10.65 ±0.34 


1.80±S;|? ±0.25 


> 3.13 


1.7+^ 


9.18 ±0.37 


0.22t^ 4 ±0.03 


Total 


86.0+-;° 




H.2 4 :^ ± 1.6 



in Fig. |^ (b) and in Table ||. Summing up the par- 
tial branching fraction for each bin, we obtain the to- 
tal B° — > D°pp branching fraction B(B° — > D°pp) = 
(I.I2+Q J4 ± 0.16) x 10~ 4 . We apply a similar procedure 
to the D°p invariant mass. The results are presented in 
Fig. I (c). In this case the total B° — » A>°pp branching 
fraction is B(B° -> D°pp) = (l.llto.u ± °- 16 ) x 1Cp4 ' 
consistent with the previous estimate. The difference 
with the result of the simultaneous fit presented in Ta- 
ble [j] is interpreted as a model-dependent error. The same 
model dependence is assumed for the B° — > D*°pp chan- 
nel. 

We examined the possibility that other B meson decay 
modes might produce backgrounds that peak in the sig- 
nal region by means of a MC sample of generic BB events 
that corresponds to about 2.6 times the data statistics. 
No peaking backgrounds were found. 

The following sources of systematic errors were found 
to be sizeable: the tracking efficiency (2% per track), 
proton/antiproton identification efficiency (3% per par- 
ticle), kaon identification efficiency (2%), 7r° efficiency 
(4%), efficiency for slow pions from D* —> Dn decays 



(8% both for n + and 7r°), branching fraction uncer- 
tainties (2% - 6%), model-dependent error (5%) and MC 
statistics (3% for B -> Dpp, 6% for B -> D*pp). The 
tracking efficiency error was estimated using rj decays to 
77 and 7r + 7r~7T°. The proton identification uncertainty 
was determined from a sample of A — > p7r~ events; the 
error in kaon selection is obtained from D* + — ► D°7r + , 
D° — > K~ir + decays. The ir° reconstruction uncertainty 
was obtained using D° decays to AT~7r + and K^tt + tt°. 
The uncertainty in the AA? signal shape parameterization 
(3%) was determined by varying the mean and width of 
the signal Gaussian within their errors. The combined 
systematic error is 14% for B° — > D°pp (B + — > D + pp) 
and 17% for B° -> D*°pp (B+ -> D*+pp). 

In summary, we report the first observation of the 
color-suppressed B° — ► D°pp and B° — > D*°pp decay 
modes. The measured branching fractions are B(B° — ► 
D°pp) = (1.18 ± 0.15 ± 0.16) x 10~ 4 and B(B° -> 
D*°pp) = (1.20±g:l| ± 0.21) x 10" 4 with 12cr and 5.6a 
statistical significance respectively. No signal is ob- 
served in the B + — > D + pp and B + — > D* + pp final 
states. The corresponding upper limits at 90% C.L. are 
B(B+ -» D+pp) < 0.15 x 10" 4 and B{B+ -> D*+pp) < 
0.15 x 10" 4 . 
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FIG. 3: (a) Dalitz plot for B° — > D°pp candidates in the B signal region. The triangles show events in the data and the small 
dots show the phase space simulation. The corresponding invariant mass spectra obtained by fitting the A_E distribution in 
each bin are shown in (b) for pp and (c) for D°p, the data indicated by points and the phase space MC by histograms. 
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